Approximate Theory of Temperature Coefficient of Resistivity of Amorphous 

Semiconductors 
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In this paper, we develop an approximate theory of the temperature coefficient of resistivity (TCR) 
and conductivity based upon the recently proposed Microscopic Response Method. By introducing 
suitable approximations for the lattice dynamics, localized and extended electronic states, we produce 
new explicit forms for the conductivity and TCR, which depend on easily accessible material parame- 
ters. The theory is in reasonable agreement with experiments on a-Si:H and a-Ge:H. A long-standing 
puzzle, a "kink" in the experimental logjg c vs. 1/T curve, is predicted by the theory and attributed 
to localized to extended transitions, which have not been properly handled in earlier theories. 

PACS numbers: 71.23.An, 71.38.Fp, 71.38.Ht. 



I. INTRODUCTION 

The temperature coefficient of resistivity (TCR) of 
an amorphous semiconductor (AS) is not only an im- 
portant quantity in transport theory, but also a critical 
parameter controlling the sensitivity of uncooled mi- 
crobolometers employed in thermal imaging "night vi- 
sion" applications^*^. 

The conventional approach to transport coefficients 
is the kinetic method (Boltzmann or master equations 
etc.). However, this is not applicable even to crys- 
talline semi-metals and semiconductors (the so-called 
Landau-Peierls criterion)^"^. Comparing to metals, the 
low carrier concentration in these materials results in a 
lower kinetic energy of carriers. Thus neither the elas- 
tic scattering by disorder, nor the inelastic scattering by 
a phonon has a well-defined transition probability per 
unit time^— . In AS, the strong electron-phonon inter- 
action of localized states requires a reorganization of the 
vibrational configuration for any transition involving lo- 
calized state(s)^. For these intrinsic multi-phonon tran- 
sitions, the energy conservation between initial and fi- 
nal electronic states (a basic condition of Fermi's golden 
rule)'^^, is violated more seriously than that for single- 
phonon emission and absorption. 

In addition, transitions between localized and ex- 
tended states (LE and EL) are not treated adequately in 
a kinetic approach. The Miller- Abrahams theory^ and 
its extensions suppose that LE and EL transitions do not 
directly contribute to conduction, and only maintain the 
distribution of carriers between localized states and ex- 
tended states in thermal equilibrium (when an external 
electric field is absent) or in the non-equilibrium sta- 
tionary state (when an external field is present). Elec- 
trical conduction is fulfilled by the transition from a lo- 
calized state to another localized state (LL) and the tran- 
sition from an extended state to another extended state 
(EE)^"^'". The theory of phonon-induced delocalization 
and the theory of transient current excited by photon 
have heuristically estimated conductivity from LE and 
El transitions. 



Rigorous expressions for the conductivity and Hall 
mobility in AS have been obtained in the microscopic 
response method (MRM)'''^'^. These expressions require 
transition amplitudes rather than transition probability 
per unit time^i. Thus the long-time limit required in a 
kinetic approach^ is avoided. To the lowest order self- 
consistent approximation, there are 29 processes con- 
tributing to conductivity and 10 processes contributing 
to Hall mobility''. For example, in a n-doped AS the con- 
ductivity from LE transitions driven solely by an exter- 
nal field is^ 
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where the real part takes the upper sign, imaginary part 
the lower sign. Qg is the physical infinitesimal volume 
element used to take spatial average. An AS can be 
viewed as uniform when we measure its properties (e.g. 
conductivity) at a linear length scale larger thani^ lOnm. 
If we take as a sphere with a radius larger than 5nm, 
then the choice of the center s of Hg inside the AS will 
not affect^ilS a^f,. Ne is the number of carriers in the con- 
duction band inside Qg, and / is the Fermi distribution 
function. The velocity matrix elements in Eq.((l]l are de- 
fined by 
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where £^ and (p^ are the eigenvalue and eigenfunction 
of localized state A. We will use letter A with or without 
a natural number subscript to denote a localized state, 
similarly Eg and Xb the eigenvalue and eigenfunc- 
tion of extended state B. Ib^a± arise from integrating out 
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the vibrational degrees of freedom, and are functions of 

external field frequency a.: TABLE I: Parameters for vibrational spectrum 



1 V-^.u^^'^*, 



B(GPa) /<(GPa) c(103m/s) jcD(A-i) p,„{g/cm^) 



lB^^^(a;) =exp{-^^coth^^(0^)2} a-SO^^ 100 52 6.21 1.44 2.33 



exp{-^(0^)^[coth— —^^cosswa -fsinsa;^]}, (4) 
z ^ z 

where o^^g = (£^ — E^)/h, coa is the frequency of 
the a* (a = 1, 2, ■ ■ ■ SA/") normal mode, JV is num- 
ber of atoms inside Qg- Denote as the shift in 
the origin of the a* mode induced by the electron- 
phonon (e-ph) interaction in a localized state^'^ A, 6'^ = 
©^{MaOJci/K)^^'^. To make the narration specific, we 
hereafter discuss conduction band transport only. For 
transport processes in the valence band, one may repeat 
the discussion mutatis mutandis. 

To calculate conductivity strictly, one needs (i) the 
eigenvalues and eigenvectors of single-electron states 
and (ii) the eigenfrequencies and eigenvectors of the 
normal modes and the electron-phonon coupling. These 
can be approximately obtained by one step of ah ini- 
tio molecular dynamics for an optimized configuration. 
Then one can compute (i) for all localized states and 

extended states; (ii) 6^ for all normal modes in each lo- 
calized states; (iii) time integrals Ib^a^ for a given oj; 
and (iv) sum over all localized states and extended states 
Y^AB- Although the result obtained in this way should 
be accurate and predictive, it is useful to develop an 
approximate theory, which also provides fimctional de- 
pendence of transport on various material parameters. 

In this paper, we will first present a tractable model 
for the conductivity and Hall mobility in AS. Then we 
will use this model to simplify the conductivity expres- 
sions obtained in the MRM for the three simplest tran- 
sitions: LL, LE and EL transitions driven solely by ex- 
ternal field, cf. Fig. 2a, 2b and 6a of 10] • The conduc- 
tivity from EE transition caused by disorder has been 
solved in the coherent potential approximatior>i2ili, ex- 
hibits weak temperature dependence, and we will not 
consider it further. 

The outline of the paper is as following. In SecHII 
we describe our approximation for the lattice vibra- 
tions and e-ph in coupling. In Sec lIII Al we first il- 
lustrate that the MRM conductivity can be put in the 
customary form of relaxation time approximation and 
of Greenwood formula. At moderately high tempera- 
ture, we invoke an asymptotic expansion to simplify the 
time integrals 1b^a±- Under the approximations intro- 
duced in Secini one can (i) obtain the velocity matrix 
elements analytically; (ii) partially carry out the two- 
fold summations over the initial and final electronic 



states. The conductivity from EL transitions is obtained 
in Sec lIII Bl The conductivity from LL transitions is cal- 
culated in Sec lIII Cl The matrix elements of electronic 
velocity could be carried out in a spherical coordinate 
system analytically. The conductivity from the LE tran- 
sitions is the same order of magnitude as those from the 
LL transitions. Below a crossover temperature T*, the 
later is larger; above T*, the former is larger. This phe- 
nomenon is the main reason for the kink in the experi- 
mental logjQ u vs. 1/T curve. As a demonstration, the 
numerical results for n-doped a-Si:H and a-Ge:H sam- 
ples are given. 

II. APPROXIMATE IMPLEMENTATION OF MRM 

A. Vibrations 

To calculate the e-ph interaction for a localized state, 
we need the transformation matrix between the atomic 
displacements and normal modes^ . Because most amor- 
phous materials are isotropiCii*^ and only acoustic modes 
are important for the e-ph interaction in one component 
semiconductors^^, one can use the acoustic dispersion 
relation for the vibrational spectrtmi: 

Lo-^ = ck, k = |k| (5) 

where CO]^ is the angular frequency for any mode char- 
acterized by wave vector k. For every k, there are one 
longitudinal and two transverse modes. We will use kr 
to label a normal mode, where t = 1, 2, 3 is the index of 
phonon branches^^. Although translational invariance 
is destroyed in AS, standing wave modes are still well- 
defined. Here, c is the average speed of soimd: 

3 2 1 

=3 = 3 + 3' 
Cf q 

where Cf and c; are the speeds of transverse and longitu- 
dinal waves which are determined by^^ the bulk mod- 
ulus B and shear modulus j-i. The cutoff wave vector 
ku = {6n^naY^^ is determined by the number density 
«a = M /V oi atoms, where V is the volume of an AS, 
M is total number of atoms^^. Ua can be inferred from 
the observed mass density p^. For a-Si and a-Ge, p^, 
B, ku and c are listed in Table U For a-Si, the De- 
bye frequency cvd is 8.91 x lO^'^Hz, not far from the ob- 
served cut-off frequency22 70meV= 1.07 X lO^^Hz. 

It is convenient to use {x3y_i)+i,X3y_i)+2/^3(;-i)+3} 
to represent the vibrational displacement vector u. = 
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{ujx, Ujy, Ujz} for the f' atom (; = 1, 2, 3 ■ ■ ■ Af). Denote 

©a (a = 1, 2, ■ ■ ■ , SAf) as the normal coordinate of the a*'' 
mode, so that the atomic displacements and the normal 
modes are related by 



X]A„,<»0^, m = 1,2, - ■ ■ ,3A/' 



(7) 



where A is the minor of the determinant \Aji — CiP-MjSj] \ 
{],! = 1,2, 3 ■ ■ ■3A/'), A is the force constant matrix^!.. 
When we use kr to label modes, Yjk ~^ Ekr- 

For a localized state, the shifts in the origins of normal 
modes caused by the e-ph interaction are the key quan- 
tities to determine the reorganization energy for transi- 
tions involving the localized stated. The shift in origin is 
determined^ by A~^, A and the e-ph coupling constant. 

and A are complicated for a system with many 
atoms. To avoid using A^^ and find a more practical A, 
we use a continuum to model the discrete random net- 
work of AS. In a continuum one can classify the atomic 
vibrations according to possible standing wave modes. 
There is no reciprocal lattice for AS. Because a contin- 
uum is isotropic and has continuous translational sym- 
metry, the wave vectors of the possible standing waves 
(k points) is imiformly distributed in the wave vector 
space (Debye sphere Sp). The M k-points inside cor- 
respond to 3J\r vibrational modes. 

The atomic displacement u at position R and time t 
satisfies the wave equation 



1 8^u(R,0 

c2 3f2 



V2u(R,f)- 



The plane wave solution of Eq.||Hll isi^ 



u(R,0 



(8) 



(9) 



kr 



where e^x is the polarization vector of mode kr. For a 



one-component systen>i^, 

ekr ■ Ckr' = ^tt'- 



(10) 



©kx and (p^^ are the amplitude and phase of mode kr, 
and are determined by the initial conditions. The in- 
verse of Eq.l(9| is 



(11) 



The normal coordinate of mode kr is ©kx^ '('^Ic+cp]^^)^ 
that 



and 



Au(R),kT ^ Skx, 



jkT,u(R) ^ ^kr- 
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TABLE II: Parameters for electronic state 
E,(eV) U(meV) nioc(A-3) Z e qTF(A-i) b 

0.121 
0.170 



a-Si 0.5[3O] 50[31J 5/10.863[32] 4 11.68 1.7 
a-Ge 0.5 51 5/11.323 4 16 1.7 



In other words, the u(R)''^ column of matrix A^^ is 
the (kr)* eigenvector belongs to the (kr)* eigen- 
value (cckx)^ = (cA:)2 of the matrix of force constants. 
Eqs.l ll2ll3l l as consequences of Eq.||8ll is contained in the 
Debye assumption ||5|. 



B. Localized states 

To obtain analytical expressions for the e-ph interac- 
tion in a localized state and the velocity matrix elements, 
we need reasonable and simple approximate wave func- 
tions for localized and extended states. We assume all 
localized states are spherically symmetric. The differ- 
ence among localized states is expressed by the localiza- 
tion length^. For a localized state A, denote R^^ as the 
position vector of the center, the normalized wave func- 
tion is 



<^^(r-Ryi) = n 



(14) 



where r and l,^ are the coordinate of electron and local- 
ization lengths^. Following Mott, is determined by 
the eigenvalue E of localized state : 



hZe^ 

4:TZ£o£ 
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where Z is the effective nuclear charge of an atom core, 
e is the static dielectric constant. Ec is the mobility edge, 
b is a dimensionless constant, b is determined by the 
shortest possible localization length fj^^^^ with E = 0. 
Realistic calculations of tail states are given in I23l - l27ll . 

The parameters^'^2 foj- electron-core interaction and 
localized state are listed in Table [III In a-Si:H and a- 
Ge:¥L^'^^, the most localized states are associated with 
dangling bonds. The localization length is one half the 
average bond length: ^j^jj^ = 2.35A / 2 and 2.45A / 2. Us- 
ing Eq.lfTSll, one has b = 0.121 and 0.170. The measured 
value of mobility edge for a-Si is rather disperse d^V^ : 
0.2-2eV: we will take^ Ec = 0.5eV. FigH] plots local- 
ization length vs. eigenenergy, we purposely left out a 
small neighborhood [Ec — U, Ec) of E, where Lf is the 
Urbach energy for band tail. When ^£ is larger than the 
linear size of a physical infinitesimal volume elementi^ 
(~100A), the corresponding localized state acts like an 
extended state for purpose of transport. 

There is a distinction between a large polaron and a 
carrier in a weakly localized state with ^ several tens 
of A. A large polaron can move freely before meeting 
a scatterer, while a localized carrier in AS is trapped in 
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the region where (p^^ has support. To make a localized 
carrier move, thermal activation involving a reorganiza- 
tion of vibrational configuration is necessary''. 




^0 0.1 0.2 0.3 0.4 0.1 0.2 0.3 0.4 



(a) (b) 

FIG. 1: Localization length as function of energy, [l(a)[ a-Si; 
m\ a-Ge. 



For many AS'^-'--, in the range of band tail, the density 
of localized states satisfies 

J[R,e,(p;E) = N{E) = !^e-{E.-E)/(i, (17) 

where U is the Urbach energy, ni^^ is the number of lo- 
calized states per unit volume. The pre-exponential fac- 
tor is determined from the requirement that the integral 
of N{E) over all localized energy spectrum should be 
nioc- In general Ec and U take different values for the va- 
lence band and the conduction band^^. Denote n as the 
carrier concentration, the Fermi energy £f of a weakly 
doped AS is: 

Ep = Ec + Uln{n/2ni,,), (18) 



Because (i) no translational invariance exists in an AS; 
and (ii) a localized electronic state is confined in some 
finite region, the spatial distribution of localized states 
needs special attention. For various macroscopic prop- 
erties, an AS can be viewed as isotropic and uniform at a 
length scale larger thanJ^ lOnm (this effectively defines 
the physical infinitesimal volume element O). Therefore 
it is convenient to describe the spatial distribution of lo- 
calized states in a spherical coordinate system. For a 
given origin and polar axis, the sum over localized states 
Ai can be changed into an integral over a combined spa- 
tial and energetic distribution of localized states: 



E 



Rc f7T fin 

R^dR / sinOde 
Jo Jo 



d<p / dEf{R,e,<p;E), 
(16) 

where R is the distance between the origin and the cen- 
ter of a localized state (pj^ , Rc is the radius of an 
AS sample, f{R,9,<p;E) is the number of localized states 
in a volume element defined by {R, R + dR), {6, 9 + d9) 
and {(p, <p + d(p) with energy (E, E + dE), i.e. position de- 
pendent density of states. Since a volume element with 
a linear size of lOnm is representative for an AS, in the 
calculation of transport coefficients, one may replace the 
volume V of the entire AS sample with the volume Q of 
a physical infinitesimal volume element. Then R^ is the 
radius of Ci. 

In a physical infinitesimal volume Q, various possible 
atomic configurations appear according to the proper 
statistical weights which would be found in a much 
larger sample. Therefore the coarse-grained average / 
of f{R,6,cp;E) over such a physical infinitesimal volume 
element is no longer position dependent: / = N(E), 
where N(E) is the usual density of states. However the 
weight factors in Eq.|[T6ll play an important role in deter- 
mining transport properties. The reason is that although 
/ is independent of (R, 6, <p), the transition amplitudes 
(velocity matrix elements) depend on the relative posi- 
tion of another localized state or on the wave vector di- 
rection of the involved extended state. 



When n < 2w;o(,, all occupied states are localized at 
T= OK. For a-Si, the conduction band energy spectr um 
((T7|l is illustrated in Fig|21 We can see from Fig l(a)| and 
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FIG. 2: Density of states of the conduction tail for n-doped a-Si 
samples: the first three vertical lines are the Fermi energy for 
n = 10^^, 10^'' and 10^^ cm^^. The rightmost vertical line is 
the mobility edge. 

Eq.tflTll that most localized states in a-Si have a localiza- 
tion length in the range 6-12A. In approximation iflTl l, 
the density of states N(E) of localized states reaches its 
maximum at E^. Therefore, the most probable localiza- 
tion length is f = cZe^ / {AneoeU). For a-Si, f = 11.75A. 
This is consistent with various experiments^—. 

Making use of relation (|T5l l, the integral over the en- 
ergy eigenvalues of localized states is converted into an 
integral over localization lengths: 



dEAN{E^ 



bZe^ 



nio 



Aneo^U Jo 



d^ bZe^ 



)■ 



(19) 

Comparing Eq. l(T6l l with the sum over states ~^ 
■Ibz ^2n)^ ^ crystal is helpful, where k is the wave vec- 
tor of a Bloch state in Brillourn zone, V is volume of the 
crystal. The matrix elements behind X^j^ may depend on 
the direction of k, d^k = k^dk sin 6]^d6i^dcl)^ takes into ac- 
count the dependence on the two wave vectors of two 
Bloch states. 
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C. Extended states 

If one imagines that an AS is obtained from deform- 
ing its reference crystal, an extended state in the AS can 
be viewed as a superposition of a principal Bloch wave 
with a given wave vector and its scattered secondary 
waves^2sM. The scattered waves are produced by scat- 
tering the principal Bloch wave with the disorder poten- 
tial (the difference between the potential energy in the 
AS and that in its reference crystal)StS,. Excepting the 
EE transitions driven by external field, we may approx- 
imate an extended state Xbi ('") by a plane wave with 
certain momentum p, and its eigenenergy is that of the 
plane wave: 



V 



-i/2gipT/r, 



(20) 



where V is the volume of AS sample, the energy zero- 
point of extended states is at the mobility edge Ec- An 
extended state in an AS is labeled by the wave vector of 
its principal Bloch wave. The sum over extended states 

becomes an integral over momentum: J^g^ J p^^-^^a ■ 



where g^^j is the e-ph coupling constant in state (pj^. Be- 
cause we consider only localized state (p^, it is conve- 
nient to shift the origin of coordinate to the center Ry\ of 
(py^. In Eq. (|22] |, we sum over all the atoms in V of an AS 
sample. In addition, the factors y(r — R„) and ^^(r — 
Rj\) in the integrand of Eq.||25]l involve two atoms, di- 
rectly integrating over coordinate is difficult (requiring 
ellipsoidal coordinate system). To obtain the coupling 
constant g^i^^, we Fourier transform dV{r — Rn)/dX„a in 
the LHS of Eq. l|23] |, first carry out the integral in coordi- 
nate r, then execute the integral over wave vector q. The 
final result is: 



16(ZeV47reoe) Xna Rn e'^^"/^ 
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]}, (24) 



D. Interaction between a carrier and an atomic core 

In a solid, the attraction to an electron from an 
atomic core may be crudely approximated by a screened 
Coulomb potentiali^ 



V{r) 



2g2 g-qrF'- 

4:TZ£oe r 



(21) 



where R„ 
„th 



|R« 



R4 1 is the distance between the 



atom to the center of localized state cp^. The 
first term decays exponentially, the second and the third 
term contain additional decay factors R^^ and re- 
spectively. Since we are concerned only with localized 
state (py^, hereafter we drop the subscript A on ^ and g. 



F. Polaron formation 



where r is the position of electron relative to an atomic 
core. qjF = 2.95{rs/a())^^^^A^^ is the Thomas-Fermi 
wave vector, Ts /uq is a number about 2 to 6. For a-Si:Hi 
and a-Gei^, we take the value for c-Si and c-Ge: qjp = 
1.7A-1. 



E. Electron-phonon coupling in a localized state 

We consider the mean e-ph interaction in a localized 
state (p^ . The e-ph interaction Hamiltonian is 



e-ph 



^ ay(r-R„) 



(22) 



where R„ (X„x, X„y, X„z) is the position vector of the 
n^'' atom, u„cr is the cr^^' Cartesian component of vibra- 
tional amplitude of the n*'' atom. Usually the average e- 
ph interaction in state (p^^ is written in a linear coupling 
forir>^ 



d^xcp\{i - RA)He_ph</'A(r - 



'■'^nagna' 

na 

(23) 



The static displacements of atoms induced by the e-ph 
interaction measure the strength of e-ph interaction and 
determine whether the e-ph coupling should be treated 
as a perturbation or be included in the zeroth order 
Hamiltonian^. The static displacement of the m^^ atomic 
degree of freedom caused by the e-ph interaction in lo- 
calized state <py^ 



.AO 



E(A- 

V 



m,p = 1,2, ■ ■ - SA/", 



(25) 



where A~^ is the inverse of force constant matrix. The 
shift 0^ in origin of the a*'' (a = 1,2, ■■ ■ 3Af) mode by 
the carrier localized in state cpj^^ 



is^ 



„A0 



(26) 



This has the physical interpretation of the polaronic re- 
laxation due to the e-ph coupling. 

If A-^andA-l were known analytically, we could use 
Eq.(|25ll to find {x^°}, and then use Eq.||26j to find {0^}. 
The continuum model in Sec lIIAl allows us to first find 
the shifts in origins {0^} of normal modes in a localized 



6 



state. Then static displacements {^m°} can be obtained 
from Eq.llTll. In the continuum model, the normal modes 
are labeled by wave vectors k. Substitute Eq.l l25|l into 
Eq.||26ll, notice A^^A-i = W^^A^, where (W^i^ 



)ixli 



Sa^Mi^ LOi^ , one concludes that 



(27) 



where n = 1,2,3- ■ - M and a = x,y, z. Substituting 
Eq.|[l2]| into Eq.l|2Zll and repl acing the sum by an integral 
over all space, Eq.l l27|l becomes 







kT 



ik-R 



(28) 



where Oa = V /J\f \s the average volume occupied by 
one atom. For a-Si and a-Gei^, Cia w (5.43A)3/4 and 
(5.66A)3/4. Eq.l[28ll expresses the shift in the ori- 
gin of normal mode k with the e-ph coupling constant 
Sncr- take k as the polar axis (z axis) and transform to 
a spherical coordinate system, because the integrand of 
Eq. l|28] | does not contain azimuthal angle <p, {g^x} 



{Sny} do not contribute to 0^.^- Only when has a 
component along k, does it contribute to 0^.^- The in- 
tegrations over the R^^ and R^^ terms in Eq.l l24|l are 
purely imaginary, and do not contribute to 0^ . The ori- 
gin shift of mode kr induced by the e-ph interaction in 
localized state is: 







kT 



1^ nZe^ / [A^neoenaf 



(2/^2] [(2/^2 + 



(29) 



Because we take AS to be an isotropic continuous 
medium, 0^.^. depends only on the magnitude k. The 
k~'^ divergence in Eq.l|29]l when k ^ is caused by 
the Debye spectrxmi (a;^ = ck). In a Debye model, 
the number of modes per unit volume per unit angu- 
lar frequency interval is^^ {2n^c)^^3k^ when k < kjj. 
The shift is smaller for higher wave number, decays 
with wave vector k as [(2/^)^ + k^]~^. Because for all 
materialsi^ qjp ~ 1.2 - 2.lA~^ while f > 2A for local- 
ized states caused by topological disorder—, the factor 
[4f - (2/^)^] in the denominator of Eqs. ll24l29l31l34l37t 
will not lead to a divergent result. 

Eq. l|29] | exhibits two obvious features: (i) 0^ > 
for every mode k; (ii) if < then 0^^ > 0^^ 
for every mode k. We have shown that three-state 
conduction processes which are first order in residual 
interactions, are the same order of magnitude as the 
two-states processes discussed hereZ.. Also, in the low- 
est order self-consistent approximation, three- and four- 
state processes must be included in the Hall mobility 
calculation'^. Some of the aforementioned transport pro- 
cesses involve at least two localized states. To carry out 
asymptotic expansion at high temperature for such pro- 
cesses, the features (i) and (ii) are essential. 
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FIG. 3: |3(a)| Static displacements x'^{t) of atoms in a localized 
state (pj^ as fimction of the distance r to the center of (p^ in a-Si: 
solid line is for f = 11.75A, circle line is for f = 23.5A. |3(b)| 
The binding energy caused by e-ph interaction as fimction of 
localization length 



The static atomic displacements in localized state A 
can be foimd from Eqs. (|26l29|l : 



/Vl/2o 3 



(27t)3 



T=l 



d^ke'^-^&^^ei,. 



(30) 



Next, substitute Eq.(|29) into Eq.||5Dll and carry out the in- 
tegral. One finds the displacement x*^^ along the radial 
direction for an atom at R caused by e-ph interaction in 
a localized state: 



R) 



4 Z*eV47Teo£ 



(31) 



where we have let /cp oo to obtain an analytic re- 
sult. It is interesting to notice that Eq. llSTI l is similar to 
the wave function of larg e polaron in strong coupling 
limit, cf. pp513-523 of ill. 

Fig 3(a)| is an illustration of Eq. (|3T] | for a-Si at ^ = 
11.75Aand 23.50A (5 and 10 times bond length). We 
observe that the more localized (smaller ^) the state, the 
larger the atomic displacements, i.e. the stronger e-ph 
interaction (larger atomic displacements). This agrees 
with previous experiments and simulation s^^ f*^^ . For 
the hardest modeSS cv = 70meV of a-Si, the amplitude 
Aq ={h/ Moj)^^^ of zero-point vibration is 0.046A, the 
amplitude Afj, ={kBT / Mco^)^^^of thermal vibration at 
300K is 0.028A. Considering these two peaks of the a- 
Si phonon spectrum are at^" 20 meV {Aq = 0.086A, 
Af,, = 0.098A) and 60 meV (Aq = O.OSOA, Af,, = 
0.033A), the static displacements of atoms estimated in 
Eq.||3T]| are twice the amplitude of vibrations. Compar- 
ing the root mean square of bond length fluctuation 0.2A 
(geometric disorder) from ab initio molecular dynamics 
simulation^'', the approximate acoustic dispersion rela- 
tion Eq.||5]l somewhat overestimates the long wave con- 
tribution in Eqs. l|28l29l30l3Hl . 
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G. Reorganization energy 

Unlike a carrier in an extended state, a carrier in a lo- 
calized state is confined by the disorder potential. Be- 
yond that, the e-ph interaction produces - an additional 
binding energy to a localized carrier in (p^ : 



A\2 



(32) 



Because the reorganization energy measures the energy 
shift from initial vibrational configuration to the final vi- 
brational configuration, is the same as^ the reorgani- 
zation energy Ab^^ of LE transition cpy^ — ?■ Xb ^^'^ 
reorganization energy A^^g of EL transition Xb ~^ 'Pa- 
AyiB = Ag^ = E^. For the continuous medium model, 
the sum over modes in Eq.(|32t may be converted to an 
integral over the Debye sphere in spherical coordinate 
system (fc, e',(p'): 



A 



AB 



2(27t)3 



dkk^x 



(33) 



where E^^ is obtained from Eq.|[34t by replacing f with 
^l, is the localization length of 



a^Oi^^ is the interference term: 



1^ nZ* / ^ne^e iJ nZ* I ^ne^t 



47r 
Mi 



^{[(2/^2)2 -(2/?i)2]-2[|i 



tan 



8(fc2,^? + 4) 16 



+ ^tan-i^ 



-2[(2/?2r-(2/^ir]-'[^tan 



tan 



Eqs. 1 1341361371 determined the reorganization energy 
A^jA] for LL transition — ?■ 



''de'sine' /'''^-'A^^2,2,^A ^2 
JO 



d<^'Mic2c2(0i^j2_ 



Owing to spherical symmetry in Eq.(|33), the direction of 
polar axis is arbitrary. Substituting Eq. (|29)l into Eq.(| 
and carrying out the integral, one finds: 



III. CONDUCTIVITY FROM LE AND EL TRANSITIONS 
DRIVEN SOLELY BY FIELD 

In this section we assemble the approximations of the 
proceeding section to estimate the various contributions 
to the conductivity. 



A 



AB 



2n- 



(2^7TZ*e2/47reo£)2g 
27Mc2Q„[(^qrf)2-4] 



(34) 



A. LE transitions driven by field 



|^tan-lM + Milium X 
Us*'''' 2 + 12 



([l + ( 



4[^ + ( 2 



15 



)}■ 



Fig |3(b)| displays the change in binding energy with lo- 
calization length. We can see that more localized states 
have larger binding energy. In other words, when a car- 
rier leaves or enters a more localized state, the required 
reorganization energy is larger, the corresponding LE 
and EL transitions are more hindered. 

The reorganization energy A^^Ai for LL transition 
(pj^^ —7- (pj^^ satisfies a reciprocity conditioi^ ^AiA2 = 
^AiAi, where 



lEM«a;2(0^-0^)2. 



Eq. (|35] | can be expressed as 



|E^I + |E^I 



M,Ai 



(35) 



(36) 
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-0 . 1 
0.15 
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FIG. 4: Conductivity and TCR as functions of temperature in 
two n-doped a-Ge:H samples at w = 0. The experimental val- 
ues are taken from jiTlliil . 
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1. Connection to relaxation time approximation and 
Kiibo-Greemvood formula 

Inside the summation of Eq.l[T]l, only electronic de- 
grees of freedom appear. Each term can be written as: 



where n = Ne / Cis is the carrier density, 

T^^(a;) =Imi[/BA+±JBA-] 
may be viewed as a relaxation time, and 



(38) 



BA\-1 



tiA\- 



2(EV 



(39) 



may be interpreted as the inverse of effective matrix ten- 
sor for transition cp^ — )• Xb- this sense, <^^jf^{(^) is 
a generalization of the energy dependent conductivitj*^ 
(r^^{cv). With this notation, Eq.((lJ becomes 

a,i,{cv) = Y: ^rfi^(a;)[l -/(EbJ]/(Ea), (40) 

a generalization of Kubo-Greenwood formula, Eq.(2.11) 
of 013 • This shows how a kinetic approach may be 
properly generalized to AS. 



2. High temperature approximation of the time integral Iba± 

To calculate Iba± (<^) defined by Eq.(|4]l, we change the 
integration variable from s to t: s = t — ifih/2. Eq.lUll 
becomes 



Iba± [c^) = exp{ - - coth \ 



^)2|g/3r,{±a;+a;^B)/2 



where 

y± n — r^r ' '^BA - l^n<^oc{aci ) ' (43) 



and 

My±) 



^n=0 n\(ln+l)' if 3/± < 1 



(44) 



The applicable condition for a-Si is T> 232K224^1<^; for 
a-Ge is T> IISK^"^. 



3. Velocity matrix elements 



Under the approximations in Sec lIIBI and III CI the 
velocity matrix elements in Eq.||2]l can be obtained by 
changing the integration variable from r to r' = r — R^, 
and introducing a spherical coordinate system with R^^i 
as the origin and p as polar axis. One can show that 

^B^A ~ ^BiA ~ ^' i-^- f^i" velocity components per- 
pendicular to p, the matrix elements are zero: 



= 0. 



(45) 



The matrix element of (the velocity component paral- 
lel to p) is 



m 



(l + p2^Vft')2 



Similarly 



w 



/ ' dfe«(±-+-AB) exp{l y{e^f csc/i^ COS fa;,}. 



w 



AB, 



-87T 



m 



2;r2 /*,2^2 



(1 + p2^Vft 



(46) 



(47) 



(41) 

If we view f as a complex variable, the saddle point 

of J Efl:(^ff )^ csc h^^ COS tcvg is at (0,0). Because the 
integrand in Eq. ll41l l is analytic in the whole complex- 
t plane, we can deform the integral path from (—00 + 
i^h/2,0 + i^h/2] to a new path C1-I-C2+C3 crossing the 
saddle point (0,0), where Ci: (—00 + ifih/2, —00 + iO], 
C2: (-00,0], C3: (0-FfO,0 + z^?!/2]. Because of the ex- 
ternal field and residual interactions being adiabatically 
introduced^, the integration along Ci is zero. When 
^bT > (oJ is the frequency of the first peak in 
phonon spectrum), YLai^aY^B^ /Kco is large. The inte- 
grals along C2 and C3 can be asymptotically calculated 
by the Laplace method^*^. The final result for 1b a± is 



Substitute Eqs. ll46l471 into Eq.lgSll, and the inverse of the 
effective mass tensor becomes 



/ BiA-1n 



m V«P/3 32nfV-'^ 



'^ff m(E0^-£0)(l + p2^V?z2)4- 

Since for each Cartesian component^, 

yt'A - cb) 



(48) 



from ||45] | and ll46l l, one has 

^Xb\^i\<PaI =Oand {^XbV\\<Pa) 



(49) 



(50) 



'BA± 



(w) = z'S/Aba 



(42a) 



(a) 



(b) 



(a) 



(b) 



o 
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FIG. 5: Conductivity and TCR as functions of temperature in 
two n-doped a-Ge:H samples at a; = 10^'' Hz. The experimen- 
tal values are taken from li^li^l . 



FIG. 6: Conductivity and TCR as functions of temperature in 
two n-doped a-Si:H samples at w = 10^^ Hz. The experimen- 
tal values are taken from 1531 



4. Relation to kinetic method 

Because (p^ vanishes at = ±oo (a = x,y,z), by 
means of partial integration, one can show that ry^g = 
— I'g^. Then 

(51) 

the last step is correct only for a cubic or isotropic body. 
For such a body, the product of two matrix elements is 
a real number. From the requirement that Re and 
Im (7a|3 are real numbers, we only require 

Re[lBA+ ± Iba-] = .y, (52) 

[e ^BA J dze iW^^ ^BA J ] 

in expression The temperature dependence | |52|| is 
the same as that obtained from the kinetic method^, 
although the Landau-Peierls condition is not satisfied. 
This is a coincidence caused by two factors. First, for 
LE, EL, LL and EE transitions driven by external field, 
the contribution to conductivity has the form of Eq.(|l]l. 
Thus only the real part of the one dimensional time inte- 
gral plays a role. In contrast to Eq.©, in the correspond- 
ing kinetic expression;^, the upper limit of time integral 
is 00 (long time limit) rather than 0. Second, because 
in both cases we apply an asymptotic expansion to cal- 
culate the time integral at high temperature, at leading 
order, the real part of (SJl is half the corresponding time 
integral in kinetic theory. The difference in temperature 
dependence only appears in subdominant terms. 



When transfer integrals or e-ph interaction are in- 
volved at first order, various transport processes are 
the same order of magnitude as the processes discussed 
here (zero-order in residual interaction). In these first- 
order processes, it is the imaginary part of a two-fold 
time integral that contribute to conductivity cf. 0]- 
Some of these first order processes do not appear in ki- 
netic models. Even for the processes expected from ki- 
netic theory, the temperature dependence derived in the 
MRM is different from that derived from kinetic theory. 



5. Summation over electronic states 

To carry out the sum over the final extended states 
and average over initial localized states, we first carry 
out X^B for a fixed localized state A. We take the center 
Ry^i of ^^(r — Ra) 3S the origin of coordinates, the in- 
cident direction k/|k| of electromagnetic wave as polar 
axis (z axis), the directions (£1,62) of two linear polar- 
ization vectors as x and y axis respectively The incident 
field is expressed as 

F = Fiei +F2e2 + Ok/|k|. (53) 

Consider an extended state (a wave packet propagating 
along p) y^i/^gip r/fi^ here for simplicity we neglected 
other waves with wave vectors close to p. We can select 
an orthogonal frame (1, m, n), where n = p/|p|, 1 and 
m are two unit vectors perpendicular to each other and 
perpendicular to n. The position vector r of electron can 
be resolved as 

r = rj_il + r_L2m + r|| n. (54) 
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According to Eq. l|50] l, one has 

(XbM'Pa) = ^{Xb\^\\\<Pa)- 



The matrix elements of the perturbation of external field 
is simplified to 



{Xb\^-^\4>a) = F-n(^B|r|||<^^). 



B. EL transitions driven by external field 

(55) Since the field-matter coupling is Hermitian, the cor- 
responding expressions for EL transition driven by field 
can be obtained from those for LE transitions driven by 
field through exchanging the status of (pj^ and Xb- 

(56) f Re Arzf bZe^tiig, 8ne^ 



Im '^^^^ 3 4n£oeU Snh^m^ Jo 



/ dp 

Jo 



sin (f 1 cos (p + F2 sin (p) - 



ihv 



BA 



where 6 is the inclination angle of p relative to k, ip is 
the azimuth angle of the orthogonal projection of p on 
plane (ej, £2) relative to ei. In this coordinate system. 



V 



(27Tfi)3 Jo 



dpp^ 



dOsine I d(p. 



(57) 



The incident field | |53] l has only x and y components. 
So that only the xx, xy, yx and yy components of the 
conductivity tensor are involved in the conduction pro- 
cess driven by field | |53l l. In consonance with Eq.|[56ll, one 
should make the substitution 



[1-/(E^)]/(Eb)- 



p4 ^exp( 47^f(jj, 



(60) 



2(/cbTA^b)1/2 
where 



\e ^ 



For the LE transition driven by the transfer integral and 
the EL transition driven by e-ph interaction, one does 
not have this symmetry^. 



— > v^^^ sin d cos (p, ~^ 



'BA 



BA 



'W 



sin 9 sin (p, (58) 



in the conductivity tensor (|T]|. The angular part of inte- 
gral l|57t can be carried out. From Eqs. 1 157158b , one can 
see axy = (Jyx = and Cxx = o'xjy = c^- Because the fac- 
tors in Eq.lUI do not depend on the position of localized 
state (pj^, one can carry out the spatial integral in 
The conductivity from LE transitions is 



Re 
Im 



4/1^ bZe^njgi. 8ne^ 
3 ineoeU Snh^m'^ Jo 



d^ / dp 

Jo 



[1-/(EbJ]/(E^)- 



?exp(- 



_bZiP_\ 



(l + p2^Vft')4 



(59) 



2(fcBTABA)i/2 



BA (2 I '''"BA+'"'" -)2 



where n = Ne/Qg is the carrier concentration, Ea and 
Ebj are given in Eqs.| |15|20] |. From Eq.|[59]l, one can eas- 
ily compute TCR: p^^^ 



dT 



-a an important ma- 



terial parameter for bolometeri«SS,. a and TCR are ex- 
pressed with easy access quantities: U and Ec for local- 
ized states, e and qjF for the interaction between elec- 
tron and atomic core, the averaged sound speed c for 
the vibrations. 



C. LL transition driven by external field 

One can similarly find the conductivity from the LL 
transitions driven by external field (Fig.2a of 01): 



Re 
Im 



2n7 



^Im 



AA, 



{^'^\A,-^\A){<,Ay 



i[lA,A+ ± Ia,a-] [1 - /(E^ J]/(E^), (62) 
where the velocity matrix elements are 

""^^1 = "!^/ ^^^'^('■-RA)3|^'/'*(r-RAi), (63) 



and 



^A,A = -? / d3xr(r-RAj^'|'(r-RA). (64) 



m 



dXa 



v'\^j^ is given in Eq. l|A2| l and w^^^ = —^a^^a- 
integral 



1 



I^^^±(a;) =exp{-^^(0^ -0^)2cothP^} (65) 



j: 
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3M2 



-(coth 



2 ^ 2 



COS (jJaS + i sin WaS)], 



contains the primary temperature dependence of con- 
ductivity. At high temperature fcgT > fiZJ, Ia^a± re- 
duces to 



-mix 



(66) 



where 



and 



AiA _ [K^<^ + <^AAi)? 

y± - 



i^A^AksT 



(67) 



(68) 



To carry out the summation over initial and final elec- 
tronic states, we first fix the initial electronic state A. 
We take the center of localized state (^^(r — R^^) as 
the origin, the incident direction k of the electromag- 
netic wave as the polar axis. Denote R = RaAi = 
|R^j — R^l the distance between the centers of local- 
ized states cpj^^ (r — R^j) and cp^^, the unit vector along 
(RAi-RA)isnA/ii = (X^^j,Yyi^j,ZAAi)/-RAAi/ where 
{XAAi'^AA^fZ'AAi) are the Cartesian components of 



vector R - 



R 



A- 



r = r±ilAA, + r±2T^AAi + r\\nAAi, 



(69) 



where IaAi and ttvaai are two unit vectors perpendicu- 
lar to each other and to n^^^^. From Eq. l|A4|l , one has 

{(Pa,\^\^a) =^aa,{'Pa,\^\\\'Pa)- (70) 
By means of Eq.([70)l, the perturbation of the electric field 



IS 



(71) 



ihv„ 



: sin 0(Fi cos (P + F2 sin <^) — L 

(La — t 



where ^ has been obtained in appendix [A] The an- 
gular integrals in summation Y^Ai can be effected: <Txij = 
(Jyx = and (Txx = CTyy = CT- Because of the uniformity 
of AS, the spatial integral in YjA can be carried out. The 
conductivity from LL transition driven by field is 



Re 
Im 



3 UneosU^ Jo i\ 



— exp(- 



hZe^ 



47reoeLf?i 



§exp(- 



hZe^ 







47reoe!J^2 



)ne' 



[l-f{EA,)]f{EA) 



2(A. ^fcBr)i/2' 



2(E0 
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FIG. 7: Conduction as function of frequency at T=300K. |7(a)[ 
LL transition; |7(b) | LE transition 

Since the conductivity tensor is usually expressed in 
a system of Cartesian coordinates, we introduce an aux- 
iliary Cartesian system (ei, €2, k), where ei and 62 are 
the two linear polarization vectors. The electric field F 
has only x and y components: F = F^ei + F2e2 + Ok. 
Because we sum over A^, the centers R^^ of localized 
states (pA^i^ ~ ^Ai) sit at different points. To simplify 
the calculation of the velocity matrix elements, we re- 
solve the position vector r of electron in an orthogonal 
frame: 



3 Ml 



AiA\/ AiA\* 
^" )(^|| ) 



'II 



(72) 



where and in appendix, to shorten the symbols, we use 
^2 instead of ^Ai' instead of ^a- 

We can see from Eqs.l l59l60l72|l that when cv = 0, 
Im (7 = for LL, LE and EL transitions. For two n-doped 
a-Ge:H samples with n = 10^^ and lO^^cm"^, log^Q a 
and TCR from LL, LE and EL transitions as functions 
of temperature at a; = are plotted in FigH) The cor- 
responding results at a; = 10^'^ Hz are plotted in Fig. 
|5l Re cr increases with frequency while TCR decreases 
with frequency. For two n-doped a-Si:H samples, the 
conductivity and TCR as functions of temperature at 
CO = 10^'^Hz are plotted in Fig[6l the results at a; = was 
reported in t54j|. The calculated TCR for a-Si:H falls^^ in 
the observed^— range between -2% and -8%. 

At a; = 0, the conductivity from LE transition is the 
same order of magnitude as that from LL transitions, the 
conductivity from EL transitions is much smaller than 
those from LL and LE transitions. There is a crossover 
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temperature T*, below T* the conductivity from LL tran- 
sitions is larger than the conductivity from LE transi- 
tions, above T* the conductivity from LE transitions is 
larger. Because the activation energy for LL transitions 
is different to that for LE transitions, this phenomenon 
explained the kink on the observed log^pU vs. 1/T 
curve . 

For two n-doped a-Si:H samples at 300K, log^QO'{co) 
vs. log^^g a; in a frequency range 10^'^ to lO^^Hz is illus- 
trated in Fig. |7(b)[ We can see that (i) the conductivity 
of LL transitions slowly decreases with co; (ii) the con- 
ductivity from LE transitions increases rapidly with fre- 
quency. The total conductivity is a sum from various 
processes'', and the conductivity from LL transitions is 
smaller than that from the LE transitions at higher fre- 
quency. The total conductivity arises mainly from LE 
transitions at higher frequencies. The general trend in 
log-j^Q cricv) vs. logjQ CO is not far from Tanaka and Fan's^ 
result cr[co) ~ co'^, but obviously deviates from the sim- 
ple power law around lO^'^Hz. We must be cautious that 
the results derived in this work is only suitable to the 
contributions from electrons: at such high frequency the 
ionic contribution should also be included. 



IV. CONCLUSION 

The microscopic response method expresses transport 
coefficients with transition amplitude rather than transi- 
tion probability per unit time, and may be used in amor- 
phous semiconductors in which Landau-Peierls condi- 
tion is violated^'-. 

We presented an approximate theory for the conduc- 
tivity and Hall mobility in amorphous semiconductors 
systematically derived from the MRM. We obtained the 
temperature dependence of the conductivity from the 
three simplest transitions: LL, LE and EL transitions 
driven solely by field, cf. Eqs. l|62l59l60l l. The conduc- 
tivity is expressed in terms of accessible physical quan- 
tities: mobility edge, Urbach energy, static dielectric 
constant and elastic modulus. LE transition (ignored 
in previous theories) contributes to conductivity in the 
same order as LL and EE transitions. Below a crossover 
temperature T*, the conductivity from LL transitions 
is larger than that from LE transitions; above T, the 
conductivity form LE transitions is larger. This phe- 
nomenon, and different activation energy for LL and LE 
transitions is the reason for the kink in the observed con- 
ductivity vs. 1/T curve. We show how a kinetic theory 
of transport can be properly generalized for AS. 
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Appendix A: velocity matrix elements between two 
localized states 



To calculate the velocity matrix elements in Eq. l|64t , it 
is convenient to adopt a system of spherical coordinates. 
We take the center of localized state (p/^{r — R^) 
as the origin = 0, the connection line R^^ — R/^ 
between the centers of two localized states as the po- 
lar axis. Denote r = |r — R/i| and r2 = |r — RaJ = 
[r^ + - IrRcosO]^/^, where K = R^, - Ra, & is the 



angle between r — R^ and Ry^^ 
ment can be written as 



Ra. The Vz matrix ele- 



ih 
m 



r^dr 



rTT r27T 

J sin0d0 J d(pe 



dz 



and one has similar expressions for the matrix elements 
of Vx and Vy. Because r2 does not depend on the azimuth 
angle (p, 



^A-,A 



' A^A 



0. 



(Al) 



We condense them as v ^ = 0: the matrix element for 
any component of velocity perpendicular to the connec- 
tion line between two localized states is zero. 

The (j) integral is immediate, the remaining r and 9 in- 
tegrals in c^^^ can be calculated by changing the inte- 
gration variable 9 to r2 for a fixed r. With the help of 



cosf 



2j.i^ - and sin 9d6 = the integral over 

9 becomes an integral over r2- One first carries out the 
integral over r2, then carries out the integral over r. For 
the velocity component parallel to the connection line 
between two localized states, the matrix element is 



At A 



1-17 ^ 



m 



l;r-3/2;r-3/2 



in 

m 



R R^' ?i 



R R^' ^1 



^ + 4 ) . ^ [2 - (i^Vr + 2R/r + 2)e-«/? ] 



We thank for support from the U.S. Army Research 
Laboratory and the U. S. Army Research Office under 



■[1-{R/^" + l)e- 
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Eq. l|A2l l displays the exponential decay of velocity ma- 
trix elements with distance R between two localized 
states. In the variable range hopping argument^, only 
the exponential decay of transfer integral with R is 
treated. In a process which is first order in transfer inte- 
gral, that is not discussed here, one may expect interest- 
ing new features. 

Because for each Cartesian component. 



i^{<pA,\'"'^\4'/ 
(Ea-Ea,) 



x,y,z, 



(A3) 



where ^' and ^" are defined by 



(A2) 



from dAlb and (IA2I I, one has 

{'Pa^ kx I'/'a) = and (<|)^^ |r|| |</)^) 



ihv 



AiA 

II 



(Ea-Eai)' 

(A4) 
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